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ABSTRACT 

We report the first detection of the v — 1-0, i?(0) ro-vibrational transition of HF at 2.499385 /mi 
arising from interstellar gas. The line is seen in absorption toward 3 background sources — HD 154368, 
Elias 29, and AFGL 2136 IRS 1 — all of which have reported H 2 column densities determined from 
observations of H2. This allows for the first direct determination of the HF/H2 abundance ratio. We 
find values of N(HF)/N(H 2 ) = 1.15 x 1CT 8 and 0.69 x 10~ 8 for HD 154368 and Elias 29, respectively. 
The sight line toward AFGL 2136 IRS 1 also shows absorption from the v = 1-0, -R(l) transition of 
HF, indicating warm, dense (tih > 10 9 cm -3 ) gas, likely very close to the central protostar. Ascribing 
portions of the HF absorption to warm and cold gas, we find 7V(HF)/A r (H 2 ) = (1.7-2.9) x 10~ 8 and 
(0.33-0.58) x 10~ 8 for the two components, respectively. Except for the warm component toward 
AFGL 2136 IRS 1, all observed HF/H 2 ratios are well below N (RF) / N {R 2 ) = 3.6 x 10" 8 , the value 
predicted if nearly all gas-phase fluorine is in the form of HF. Models of fluorine chemistry that 
account for depletion onto grains are able to reproduce the results toward HD 154368, but not in the 
cold, dense gas toward AFGL 2136 IRS 1 and Elias 29. Most likely, some combination of simplifying 
assumptions made in the chemical models are responsible for these discrepancies. 
Subject headings: astrochemistry - ISM: molecules 



1. INTRODUCTION 

Fluorine is unique among all elements in that its neu- 
tral, atomic form reacts exothermically with H 2 because 
the dissociation energy of HF is larger than that of H 2 . 
The ionization potential of F is larger than that of H, 
such that atomic fluorine is expected to remain in neutral 
form throughout most of the interstellar medium (ISM) . 
In regions with H 2 , fluorine will convert to molecular 
form via the reaction 



H 2 



HF + H. 



(1) 



HF is destroyed by photodissociation and ion-neutral re- 
actions (with, e.g., C + , Si + , He+, H3 ). While these 
processes — most importantly reactions with C + and 
photodissociation — can compete with reaction ([1]) in the 
diffuse, primarily atomic ISM to keep a significant frac- 
tion of fluorine in atomic form, chemical models predict 
that in molecular gas (H 2 /H > 1) HF becomes the domi- 
nant reservoir of fluorine (jNeufeld et al.ll2005l : iZhu et al.1 
12001 . 

Interstellar HF was firs t discovered using th e Infrared 
Space Observatory (ISO; iKessler et al.l [T996) targeting 
the J = 2-1 transit ion at 121.6973 /im toward Sgr B2 
(jNeufeld et al.l 11997ft . Given the large dipole moment 
and rotational constant of HF though, only in excep- 
tional cases of radiative pumping (as is the case with 
Sgr B2) or at high densities (jih > 10 9 cm~ 3 ) is the 
J = 1 level expected to be significantly populated. In 
most regions of the ISM, nearly the entire population of 
HF is expected to be in the J = ground state, thus 
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limiting the utility of the J = 2-1 line. It wa s not un- 
til the launch of Herschel (jPilbratt et al.ll2010D that the 
J = 1-0 transition at 243.244 fim (1232.47 63 GHz) could 
(and would) be routinely observed (e.g.. INeufeld et all 
2010t [Sonnentrucker et all l2010l: IPhillips et all 120101: 



Monje et alll201lHEmprechtinger et alll2012l ). Using H 2 
column densities estimated from other tracers such as 
13 CO and CH, these studies found relative abundances of 
N(HF) /N(R 2 ) in the (1-2) x 10" 8 range in diffuse molec- 
ular clouds, with an average value of 1.4 x 10~ 8 . 

Observations of neutral atomic fluorin e have also been 
made in several diffuse cloud sight lines (jFederman et all 
120051: ISnow et al] I2007D targeting transitions at 951 A 
and 954 A using the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE). All sight lines where F is detected have 
low molecular hydrogen fractions — /h 2 < 0.3, where 
/h 2 = 2iV(H 2 )/[iV(H) + 2iV(H 2 )]— so that F is expected 
to be the dominant form of fluorine. The average frac- 
tional abundance of atomic fluorine with respect to to- 
tal hydrogen— x(F) = N(F)/N H , where N H = 2V(H) + 
2N(H 2 )— is found to be x(F) = 1.8 x 10~ 8 , about 60% 
the sol ar system abun dance of2.9x 10 -8 found in mete- 
orites (|Loddersll2003l and references therein). 

Taking x(F) = 1.8 x 10~ 8 to be the interstellar 
gas phase abundance of fluorine and assuming condi- 
tions where both fluorine and hydrogen are predom- 
inantly in molecular form results in the prediction 
N(RF)/N(}1 2 ) = 3.6 x 10" 8 . This is about twice as 
large as the average value determined from Herschel ob- 
servations in diffuse molecular clouds, and about 100 
times as lar ge as the values fo u nd in dense molec- 
ular clouds (INeufeld et all H997t IPhillips et all l2010l: 
Emprccht inger et alll2012ll The severe underabundance 
of gas phase HF in dense clouds is thought to be the result 
of nearly a ll HF freezing out onto grains in such cold env i- 
ronments (jNeufeld et al. 2005; N eufeld fc Wolfirel l2009'). 
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and it is possible that a similar mechanism, albeit less 
efficient, operates in diffuse clouds as well, resulting in 
the factor of 2 underabundance in those regions. Still, 
all current estimates of the HF/H2 ratio rely on abun- 
dance relationships between H2 and some other molecule 
(e.g., CH, CO), and observations of that proxy molecule. 

In this paper, we present the first direct determina- 
tion of the HF/H2 ratio by observing the v = 1-0, R(0) 
ro-vibrational transition of HF in sight lines where the 
H2 column densitiy is known either from observations 
of electronic transitions in the ult raviolet (jSavage et al.l 
\WI% iRachford etaTl 12001 12001 or of ro- vibrational 
quad r upole transiti ons in the infrared (e.g., lLacv et all 
11994 iKulesal [20021. While higher lying ro-vibrational 
transitions of H F have been observed in atmospheres of 
cool stars (e.g.. iJorissen et al.lll992l iCunha et al.l 120031 : 
Uttenthal er et al.ll2008| ). this is the first detection of in- 
terstellar HF via ro-vibrational spectroscopy. 

2. OBSERVATIONS & DATA REDUCTION 

All observations were carried out using the Cryogenic 
High-resolution In frared Echelle Spectrograph (CRIRES; 
iKaufl et al.ll2004D on UT1 at the Very Large Telescope 
(VLT). Observations were performed in service mode, 
and CRIRES was used with its 0'.'2 slit to provide a re- 
solving power (resolution) of about 100,000 (3 km s" 1 ), 
and a reference wavelength of 2502.8 nm to position 
the v = 1-0, R(0) (A = 2.499385 fim) and R(l) (A = 
2.475876 //m) transitions of HF on detectors 3 and 1, re- 
spectively. The adaptive optics system was utilized in all 
cases (excluding Elias 29 and AFGL 2136 IRS 1 where 
no bright natural guide stars are available) to maximize 
starlight passing through the narrow slit. Spectra were 
obtained in an ABBA pattern with 10" between the two 
nod positions and ±3" jitter width. Details regarding in- 
tegration times and dates of observations are presented 
in Table Q] Our observing program was designed with 
the intent of removing atmospheric absor ption features 
via d ivision by model atmospheric spectra (jSeifahrt et al.l 
2010). As such, no telluric standard stars were observed 
(with the exception of HR 6508, paired with the final 
observation of HD 154368 on 2012 Aug 02). 

Raw data images were processed using the CRIRES 
pipeline version 2.2.1. Calibration techniques applied to 
science images include subtraction of dark frames, di- 
vision by flat fields, interpolation over bad pixels, and 
correction for detector non-linearity effects. Consecu- 
tive A and B nod position images were subtracted from 
each other to remove sky emission features. All of the 
images from each nod position were combined to cre- 
ate average A and B images. One-dimensional spectra 
were extracted from these images using the apall routine 
in IRAE0, and then imported to Igor Pro0. Wavelength 
calibration was performed using atmospheric absorption 
lines, and is accurate to ±1 km s _1 . Spectra from the A 
and B nod positions were averaged on a common wave- 
length scaleQ 

In order to remove atmospheric absorption features 
from science spectra, we attempted division by atmo- 

4 http://iraf.noao.edu/ 

5 http://www.wavemetrics.com/ 

6 Combination of A and B spectra is done after wavelength cal- 
ibration due to a slight (about one-half pixel) shift in wavelength 
along the dispersion direction between the two nod positions. 



spheric model spectra. However, simultaneous removal of 
both strong and weak H 2 features proved exceptionally 
difficult, such that the targeted continuum level signal- 
to-noise ratio (S/N) of 1000 could not be achieved. We 
then attempted dividing science spectra by other science 
spectra with the requirement that the expected inter- 
stellar feature of the "standard" star be shifted at least 
4 km s _1 with respect to the interstellar absorption fea- 
ture of the science target. This division procedure allows 
for interactive stretching and shifting of the standard star 
spectrum in the wavelength axis, and scaling of the stan- 
dard star intensity acco rding to Beer's law, and is de- 
scribed in iMcCalll (|2001[ ) . Removal of both weak atmo- 
spheric lines and the wings of strong atmospheric lines us- 
ing this method was found to be satisfactory for observa- 
tions toward X Oph, HD 149404, HD 152236, HD 154368 
(on all nights), and Elias 29. For AFGL 2136 IRS 1 divi- 
sion by a model atmospheric spectrum produced better 
results. For £ Oph the interstellar HF feature is blended 
with atmospheric HF absorption, and for HD 179406 the 
interstellar HF feature is blueshifted into the atmospheric 
v \ 96,3-85,4 H2O line, such that division by neither a 
standard star nor an atmospheric model is able to ad- 
equately remove atmospheric features. As a result, we 
omit these sight lines from our analysis. After division by 
a standard star or model, any remaining large scale con- 
tinuum fluctuations were removed by a custom-written 
procedure^ 

Spectra covering the i?(0) transition resulting from our 
reduction procedures are shown in Figure [TJ \ Oph, 
HD 149404, and HD 152236 show no evidence of in- 
terstellar HF absorption, while HF is clearly detected 
toward Elias 29 and AFGL 2136 IRS fl Spectra of 
HD 154368 from all three nights when observations were 
taken show marginal absorption features at the positions 
expected given the Earth's motion and known interstellar 
gas velocities. Because absorption is seen at the expected 
wavelength on all three nights, we are more confident in 
considering this a detection. Figure [2] shows spectra of 
the two dense cloud sight lines covering the R(l) transi- 
tion. Although no R(l) feature is seen toward Elias 29, 
HF absorption out of the J = 1 level is detected toward 
AFGL 2136 IRS 10. Figure [3] shows spectra in velocity 
space covering both the R(0) and R(l) transitions in all 
three sight lines where the R(0) line is detected. Inter- 
estingly, the R(0) and i?(l) line widths and line profiles 
differ toward AFGL 2136 IRS 1, possibly indicating dif- 
ferent physical conditions (e.g., temperature, density) in 
the gas from which the absorption arises, as will be dis- 
cussed below. 

3. ANALYSIS 

HF absorption feature s were fit with Ga u ssian func- 
tions as described in llndriolo fc McCall (|2012t l for 
the purpose of determining equivalent widths, velocity 
FWHM, and interstellar gas velocities. These parameters 
are reported in Table [5J along with HF column densities 
derived assuming optically thin absorption. In the case 
of HD 154368, spectra from all three nights were shifted 

7 http:/ /fermi. uchicago.edu/freeware/IgorPlugins/ 

8 The HF R(0) line toward AFGL 2136 IRS 1 is partially blended 
with a water feature that is also astrophysical in origin. 

9 Again, partially blended with an astrophysical water line. 
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to the LSR frame and averaged prior to this fitting pro- 
cedure. For AFGL 2f36 IRS I the HF R(Q) and H 2 
v\ 96,3-85,4 lines were fit simultaneously, as were the HF 
R(l) andH 2 v\ 107,4-96,3 lines. For non-detections, un- 
certainties in the equivalent width were derived from the 
standard deviation on the continuum where the absorp- 
tion lines were expected, the wavelength step between 
pixels, and the number of pixels expected to contain ab- 
sorption assuming a line with FWHM=3 km s _1 (for 
Elias 29 and HD 154368 the measured FWHM of the 
i?(0) lines were used in determining uncertainties in the 
equivalent width of the R(l) line). These la uncertain- 
ties in equivalent width and column density are also re- 
ported in Tabled 

Upper limits on AT(HF) in a single cloud component 
are taken to be 3 times the uncertainty of the HF column 
densities for non-detections. In sight lines where multiple 
velocity components are observed in CH (2 components 
toward HD 152236 and 3 components toward HD 149404; 
IGredel et al.l |l993) , total line-of-sight upper limits on 
7V(HF) are taken to be the upper limit determined for a 
single cloud component multiplied by the number of ob- 
served components. These values are reported in Table 
[H along with H2 and CH column densities from the liter- 
ature for sight lines we observed in HF, and the relative 
abundance of HF with respect to these species. Figure 2] 
shows iV(HF) versus A(H 2 ) (left panel) for the six sight 
lines in Table H and iV(HF) versus N(CR) (right panel) 
for sight lines in this study and reported in the litera- 
ture. Also shown in Figure |4] are the expected abun- 
dance relationships— N(RF)/N(R 2 ) = 3.6 x 10~ 8 and 
N(HF)/N(CH) — 1 — assuming conditions where nearly 
all hydrogen is H2, nearly all gas- phase fluorine is HF , 
and JV(CH)/AT(H 2 ) = 3.5 x 10" 8 (jSheffer et all l2008h . 
as well as abundances p r edicte d by chemical models pre- 
sented in iNeufeld et all ()2005D . 

4. DISCUSSION 

In the left panel of Figure 0] it is clear that our mea- 
sured HF abundances do not match the expected abun- 
dances for fully molecular H 2 and HF. Aside from the 
upper limits toward HD 149404 and HD 152236, which 
are close to the predicted ratio, all measured HF abun- 
dances are lower than expected. The right-hand panel 
of Figure 2] shows that the measured HF abundances 
also do not match the expected 1:1 ratio with CH, 
again with observed values below predicted values. Mea- 
sured CH abundance are consistent with the relationship 
iV(CH )/Af(H 2 ) = 3.5±l\ x 10~ 8 found bv ISheffer etHl 
(2008), suggesting that we are seeing low abundances 
of HF, rather than high abundances of CH. Our find- 
ings agree with HF/CH ratio s determined from Hersche l 
observations of both species dSonnentrucker et al J 120 lot 
iGerin et alJ2010HEmprechtinger et al.ll2012h . As such, it 
seems likely that the fully molecular fluorine assumption 
is incorrect, and that one or more mechanisms must de- 
crease the formation rate of HF, increase the destruction 
r ate of HF, or r e move H F from the gas phase. 

INeufeld et al.l (|2005l ) discussed interstellar fluorine 
chemistry in detail, and presented a simple model for a 
plane-parallel geometry cloud illuminated on both sides 
by the interstellar radiation field. In addition to study- 
ing the pure gas-phase reaction network, they also in- 
vestigated the effects of F and HF depletion onto grains. 



Select results from their study are plotted in the left-hand 
panel of Figure |4j The red dotted curve is for a model 
cloud with F depletion independent of Ay, Xuv = 10 
and 7Jh = 10 2 cm -3 . While the model is consistent 
with our upper limits, it cannot reproduce the observed 
HF/H 2 ratios in sight lines where HF is detected. In- 
creasing the density does not improve the agreement, 
as model results shift toward the HF/H 2 = 3.6 x 10~ 8 
fully molecular scenario (see Figure 7 in INeufeld et al.l 
120051 ). Similarly, increasing the impinging UV field only 
serves to worsen the agreement with our observations. 
The purely gas-phase chemical models still predict that 
most fluorine is in the form of HF over a wide range of 
input conditions, suggesting that formation and/or de- 
struction rates of HF do not vary enough to account for 
the observed abundances. 

Model results with variable fluorine depletion are also 
shown in Figure |4j and in some cases do a better job 
of reproducing observations. Blue dashed curves show 
model results for F depletion variable with Ay, xuv = 1, 
and n H = 10 2 , 10 2 5 , 10 3 , 10 3 5 , and 10 4 cm" 3 (from 
top to bottom). iSonnentrucker et al.l (|2007l ) estimated a 
density of jih = 240 cm -3 from C 2 observations toward 
HD 154368, and the model curve with nn = 10 2 5 cm~ 3 
is in excellent agreement with our observed HF column 
density. However, models with depletion are unable to 
reproduce our findings toward the dense cloud sight lines 
Elias 29 and AFGL 2136 IRS 1, drastically underpredict- 
ing observed HF column densities even with an increased 
UV fieldJ3 To determine why this is so, we first explore 
properties of the Elias 29 and AFGL 2136 IRS 1 sight 
lines in detail. 

4.1. Elias 29 

The line of sight toward the low-mass protostar 
Elias 29 passes through several different environments. 
Multiple velocity components are traced by emission in 
CO, 13 CO, C 1 8 Q, CS, C 34 S HCO +, H 13 CO+, H 2 CO, 
and CH3OH (jBoogert et all [20021) . Two foreground 
clouds at ulsr = 2.7 km s _1 and 3.8 km s _1 are 
thought to be unassociated with the protostar, while 
some combination of the molecular ridge in which the 
protostar is embedded and the protostellar envelope it- 
self is thought to give rise to emission and self- absorption 
seen at ^lsr = 5 km s _1 . Densities and kinetic tempera- 
tures of riH ~ 10 4 cm -3 and T ~ 20 K in the foreground 
clouds and fin ~ 4x 10 5 cm* 3 and T ~ 15 K in the molec- 
ular ridge ha ve been estimated fro m observations of C 18 
and HCO+ (jBoogert et al.ll200l . Warm (T - 500 K) 
CO and H 2 are seen via ro-vibrational transitions in 
absorption (at low spectral resolution, thus lacking veloc- 
ity i nformation), likely ar ising from a hot core or flared 
disk (jBoogert et alj|200(l . This warm gas is also traced 
by emission arising from vibrationally excited states of 
12 CO and 13 CO seen in high spectral re solution obser- 
vations near 4.6 fim ([Herczeg et al.l[2~011f ). Additionally, 
absorption in 12 CO from to -90 km s _1 (with respect to 
systemic) is thought to trace a wind/outflow component. 

The total H 2 column density in all three cold 

10 UV flux is giv en in terms of the "standard" UV background 
fromQDrainc (1978!), denoted as xuv = 1- 

11 The green dash-dot curve is for a model with F depiction 
variable with Ay, xuv = 1000, and rtjj = 10 4 cm -3 . 
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components — inferred from C 18 J = 2-1 emission — 
is N(H.2) — 2.9 x fO 22 cm~ 2 , in good agreement with 
(2.6±0.5) x 10 22 cm -2 determined directly from observa - 
tions of the v = 1-0, S(0) transition of H 2 (|Kulesall2002T) . 
This suggests that the H2 column density listed in Table 
[3] and used in Figure [4] is the total line-of-sight H2 col- 
umn density, and does not discern between the different 
cloud components. However, the observed HF absorp- 
tion feature is relatively narrow (FWHM=4.0 km s _1 ) 
and has a velocity centroid that favors association with 
one or both of the foreground clouds. Given the small 
spread in velocity along this sight line though, it is diffi- 
cult to ascertain which cloud(s) may give rise to HF ab- 
sorption. It is quite possible that we have underestimated 
the HF/H2 ratio by including H2 that is not associated 
with the component where HF is observed (and conse- 
quently overestimated the ratio in components where HF 
absorption does not arise). Making this assumption, only 
in the case where all HF is assoc iated with the 2.7 km s _1 
cloud (JV(H 2 ) - 5 x 10 21 cm- 2 ; lBoogert et al.ll2002l ) can 
the models reproduce iV(HF) and N(H.2), and even then, 
it is only models with constant depletion that are capable 
of doing so. 

4.2. AFGL 2136 IRS 1 

AFGL 2136 IRS 1 is an embedded massive protostar 
that shows a complex 3-lobed morphology at 2.2 /xm — 
thought to trace the cavity walls of a bipolar outflow — 
and a pol arization signa t ure indicative of a nearly edge - 
on disk (jKastner et all 119921 ; iMurakawa et al.l 120081) . 
Emission maps in 12 CO J = 1-0 and J = 2-1 and 13 CO 
J = 1-0 and J = 2-1 show massive molecular outflows 
that support this p icture, and give a sy stemic velocity of 
v lsr = 22 km s _1 (|Kastner et al.lll994l ). Absorption due 
to the v = 1-0 and v = 2-1 bands of 12 CO and the v = 1- 
band of 13 CO reveal both cold (T ~ 20 K) and warm 
T ~ 580 K) compon ents at ulsr = 26.5 ± 2.8 km s" 1 
Mitchell et aHll990l ). Warm (T ~ 500 K) gas is also 
observed in absorption from the vi ro- vi brational band 
of H2O (jBoonman fc vanTT ishocck 2003), and our own 
observations targeting HF have revealed absorption in 
the v\ and v% bands of H2O arising from more than 30 
different rotational levels with relative populations con- 
sistent with a temperature of 500 ± 20 K (Indriolo 2013, 
in preparation). 

As mentioned above, a very high density (nn ^ 10 9 ) 
and/or exceptional radiative pumping is necessary to 
maintain substantial population in the J = 1 level 
of HF, making our detection of the R(l) line of HF 
toward AFGL 2136 IRS 1 rather s urpris ing. How- 
ever, studies of CO (Mitchell et all 119901 ) and H 2 
(jBoonman fc van Dishoeckl 120031 ) suggest that the ab- 
sorbing gas toward AFGL 2136 IRS 1 is dominated by 
warm (T ~ 500 K), dense (ran > 10 10 cm -3 ) gas in close 
proximity to the central protostar, perhaps located in a 
protostellar disk. In this case, relative populations in 
rotationally excited levels of HF are expected to be ther- 
malized, and the -R(l) line is detectable. The profile of 
the HF R(l) line is in good agreement with those of most 
excited water lines (Indriolo et al. 2013, in preparation), 
suggesting that the two species reside in similar gas con- 
ditions. The HF i?(0) line, however, is narrower than 
the R(l) line and most water lines, implying a different 



origin. 

In thermal equilibrium at T = 500 K the relative pop- 
ulation in the J = 1 and J = levels is expected to 
be N(1)/N(0) — 2.7, much larger than the observed 
line-of-sight ratio of 1.03. However, this can be rec- 
onciled if we ascribe a portion of the R(0) absorption 
to cold, moderately dense gas in the protostellar enve- 
lope. Assuming that HF is in thermal equilibrium in 
the warm gas results in N(0) warm = 0.92 x 10 14 cm~ 2 
and N(0) co id = 1.50 x 10 14 cm~ 2 . Under this assump- 
tion levels with J > 2 will also be populated, such that 
N(HF) warm /N(l) — 6.5, and the total HF column den- 
sity in warm gas is N(HF) warm — 1.62 x 10 15 cm~ 2 . As 
all HF in cold gas is still expected to be in the J = 
level, iV(HF) ooW = 1.50 x 10 14 cm" 2 . 

Determining HF/H2 ratios then requires finding the 
H2 columns in both the warm and cold components. 
To do so, we set the temperature in the warm gas to 
500 K, and allow the temperature in the cold gas to vary 
from 20-80 KJ3 calculating the fractional population of 
H2 expected in the J = level for all cases. The sum 
iV(H 2 ; J = 0) warm + N(R 2 ; J = 0)coid is constrained to 
equal 3 x 10 22 cm~ 2 (total line-of-s ight column density 
in the J = level of H 2 reported by Kulesa 2002) , while 
the ratio of warm H2 t o cold H2 is assumed to be the 
same as that found by iMitchell etaLl (|1990t ) for warm 
to cold 13 CO: 2.14. Given the 20-80 K range assumed 
for the temperature of the cold gas, the total H2 column 
density in the warm component is (5.5-9.8) x 10 22 cm -2 , 
and the total H2 column density in the cold component 
is (2.6-4.6) x 10 22 cm -2 . Taking the HF column densities 
from above, we find HF/H 2 ratios of (1.7-2.9) x 10~ 8 for 
the warm gas, and (0.33-0.58) x 10~ 8 for the cold gasP 3 ! 

The HF and H2 column densities found from this anal- 
ysis are shown in Figure |4] Abundances in the warm gas 
component are nearly consistent with the simple scenario 
where all fluorine is in the form of HF. Given that H2O 
is observed to be predominantly in the gas phase in this 
region, it is likely that depletion of fluorine onto grains is 
negligible (due to the grains being warm), and the sim- 
ple chemistry mentioned in Section [1] applies. Relative 
abundances of HF and H2 in the cold gas component, 
however, remain inconsistent with any of the models dis- 
cussed herein. 

4.3. Discrepancy between models and observations 

The chemical models of Neu feld et al.l ([2005D that in- 
voke fluorine depletion varying with visual extinction as- 
sume the rate at which F and HF stick to grains is pro- 
portional to the dust density (which is proportional to 
the density of H nuclei assuming a fixed gas-to-dust ra- 
tio). HF is only removed from grains via UV photodes- 
orption, a process inhibited in cloud interiors where UV 
photons are only emitted following cosmic-ray induced 
excitation of H2 Lyman and Werner bands. As a result, 
the models predict that all HF inside of dense clouds 
should be frozen onto grains. However, dust grains in the 
inner portions of a protostellar envelope are expected to 

12 The 20 K limi t is based on the cold component found by 
IMitchell etTaLl 1T990D in 13 CO and the 80 K limit on an H 2 ortho- 
to-para ratio of 1. 

13 Note that a larger HF/H2 ratio corresponds to a lower tem- 
perature in the cold component. 
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be warmed by the central radiation source — an effect not 
included in the depletion model — leading to the evapora- 
tion of icy grain mantles. We have already alluded to this 
likely being the scenario for the warm gas component to- 
ward AFGL 2136 IRS 1. It is possible that some portion 
of the HF absorption seen toward Elias 29 also arises from 
warm gas close to the protostar, but the lack of -R(l) ab- 
sorption means we cannot attempt the same warm/cold 
decomposition performed for AFGL 2136 IRS 1 above. 
However, this mechanism still cannot explain the larger 
than expected HF/H2 ratio in the cold component of 
AFGL 2136 IRS 1 unless the grain temperature is signif- 
icantly warmer than the gas temperature. 

It is possible that some amount of the differences be- 
tween observed and model abundances stems from geo- 
metric effects. Model results presented in Figure@]are for 
plane-parallel clouds viewed face-on. For inclined view- 
ing angles, sight lines sample more material near the edge 
of the cloud, where, in the case of fluorine depletion in- 
creasing with Ay, HF has not yet frozen onto grains. As 
a result, the line-of-sight HF/H2 ratio increases for in- 
clined viewing angles. A similar effect occurs for the same 
reason if we consider multiple clouds along the line-of- 
sight. Still, these geometric effects alone cannot bring the 
model curve for a dense cloud (xuv = 1, tin = 10 4 cm -3 ; 
bottom dashed blue curve in the left panel of Figure 2]) 
into agreement with observations. 

The constant densities assumed in these cloud models 
may also contribute to discrepancies with observations. 
For the nn = 10 4 cm -3 model, HF is efficiently removed 
from the gas phase even at the edge of the cloud because 
the sticking coefficient is proportional to the gas density. 
Dense clouds are expected to be surrounded by diffuse 
material though, where HF will remain in the gas phase. 
A line of sight probing a dense cloud with a diffuse enve- 



lope would contain an HF abundance dependent almost 
entirely on the size of the envelope, because little or no 
fluorine depletion occurs there. In this way, a cloud that 
consists predominantly of dense gas can still show an ap- 
preciable HF column density if it is surrounded by diffuse 
material. 

5. SUMMARY 

For the first time, we have detected the v = 1-0, R(0) 
transition of HF in interstellar gas. Additionally, we have 
detected both the R(0) and i?(l) transitions in the inner 
envelope or protostellar disk around AFGL 2136 IRS 1. 
This warm gas associated with a massive protostar has 
an HF/H2 ratio similar to that predicted by the sce- 
nario where all hydrogen is in molcular form and all 
fluorine is in HF, while observed HF abundances in in- 
terstellar gas are much lower. Chemical models assum- 
ing increasing depletion of HF onto grains with increas- 
ing Ay can reproduce the HF/H 2 ratio observed toward 
HD 154368, but not in the cold, dense component to- 
ward AFGL 2136 IRS 1 and Elias 29. It is possible 
that some assumptions used in modeling the fluorine 
chemistry — constant densities, plane-parallel geometry, 
the lack of HF evaporation at high temperatures — all 
contribute to discrepancies between observed and mod- 
eled abundances, but it is difficult to ascertain which may 
have the greatest effect. More detailed chemical models 
and an expanded sample of HF observations may help to 
elucidate the most important factors. 
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Observed Wavelength (fim) 

Fig. 1. — Spectra toward target sight lines covering the v=l-0, R(0) transition of HF. Top panel: Observation of f Oph, used to 
show atmospheric transmission spectrum in the region of interest. The two strongest atmospheric lines are due to the v\ 96,3—85,4 and 
v l 96, 4^85,3 transitions of H2O at 2.499199 /im and 2.499706 fim, respectively. HF is clearly present in the Earth's atmopshere, and the rest 
wavelength of the -R(O) transition at 2.499385 [ira is marked by the vertical dashed line. Bottom panel: Spectra toward science targets. 
Gaps in spectra are regions where division by a standard star was unable to adequately remove atmospheric absorption lines. Solid vertical 
lines above spectra mark the expected position of HF absorption given the Earth's motion a nd previously repor ted interstellar gas velocities 
(note that HD 149404 and HD 152236 have multiple velocity components reported in CH; Grcdcl ct al. 1993). Interstellar HF is detected 
toward Elias 29 and AFGL 2136 IRS 1, and marginally detected toward HD 154368. The absorption feature in the AFGL 2136 IRS 1 
spectrum centered at about 2.4993 £im is from the v\ 96,3-85,4 transition of H2O and likely arises in gas very close to the protostar (Indriolo 
et al. 2013, in preparation). 




2.4756 2.4758 2.4760 2.4762 2.4764 
Observed Wavelength (|im) 

Fig. 2. — Spectra toward target dense cloud sight lines covering the v=l— 0, R(l) transition of HF. Top panel: Observation of £ Oph, used 
to show atmospheric transmission spectrum in the region of interest. The rest wavelength of the R(l) transition at 2.475876 fim is marked 
by the vertical dashed line. Bottom panel: Spectra toward AFGL 2136 IRS 1 and Elias 29. Solid vertical lines above spectra mark the 
expected position of HF absorption given the Earth's motion and previously reported interstellar gas velocities. The -R(l) line is detected 
toward AFGL 2136 IRS 1, but not toward Elias 29. As before, the absorption feature to the left of the HF line toward AFGL 2136 IRS 1 
is due to water near the protostar itself, in this case, the v\ 107,4—96,3 transition. 



AFGL 21 36 R(0) 




-40 -20 20 40 

LSR velocity (km s 

Fig. 3.— Spectra toward AFGL 2136 IRS 1, Elias 29, and HD 154368 showing the HF R(0) and R(l) lines in velocity space. The R(0) 
and R(l) absorption features toward AFGL 2136 IRS 1 are at roughly the same velocity, but have different profiles and line-widths (see 
Table [2J. Only the R(0) line is detected toward Elias 29 and HD 154368. Features between —15 km s _1 and —40 km s~ x are due to 
imperfect removal of an atmospheric water line. For HD 154368, these spectra are the combination of all three nights of data. Red curves 
show Gaussian fits to the HF (and H2O for AFGL 2136 IRS 1) absorption features. 
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A/(H 2 ) (cm 2 ) W(CH) (cm 2 ) 

Fig. 4. — Plots showing column densities of HF, CH, and H2. Left panel: ./V(HF) vs. N(H2) for the six sight lines reported in this study. 
Filled black diamonds mark detections of HF, and black arrows mark upper limits to 7V(HF) reported in Table[3] The red and blue s haded 
regions show a decomposition of the AFGL 2136 IRS 1 data into a warm and cold component, respectively (as described in Section |4.2| |. 
where dark lines mark the derived HF column densities, extending horizontally to show the range of inferred H2 column densities, and the 
lighter shaded regions extend vertically to mark la uncertainties in the HF column de nsities. The solid line shows HF/H2 = 3.6 X 10 -8 . 
Dashed (blue), dotted (red), and dash-dot (green) curves show results from the models of Ncufeld ct al. (2005) with different densities (njj), 
ultraviolet radiation fields (xuv)i an d fluorine depletion schemes. The red dotted line is for a model cloud with F depletion independent of 
Ay, XUV = 1) an d n-H = 10 2 cm -3 (see their Figure 7). Blue dashed curves show model results for F depletion variable with Ay, xuv = 1, 
and rin = 10 2 , 10 2 ' 5 , 10 3 , 10 3,5 , and 10 4 cm -3 from top to bottom (see their Figure 9). The green dash-dot curve is for F depletion variable 
with Ay, Xuv = 1000, and njj = 10 4 cm -3 (also their Figure 9). Right panel: 7V(HF) vs. iV(CH) for sight lines reported in this study, 
as well as sight lines where both HF and CH have been observed with Herschel. The filled black d iam onds and arrows as are before, while 
the open black sq uares mark HF and CH colu mn densities reported in |Sonncntruckcr ct al. (2010) and Gcrin et al.l pOlOh . and open black 
triangles are from Emprcchtingcr ct al. ( 20 121) . Th e solid line shows the ~ 1:1 ratio of HF to CH expected given fully molecular fluorine 
and the CH/H 2 ratio found bv lSheffer et alj 12008ft . Elias 29 and AFGL 2136 IRS 1 are not in this panel as no CH column densities are 
reported in the literature. 
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TABLE 1 Log of Observations 



Target 


Date 


Integration Time 

GO 


HD 154368 


2012 Jun 26 


1200 


HD 154368 


2012 Jul 04 


1200 


HD 154368 


2012 Aug 02 


1440 


X Oph 


2012 Jul 04 


240 


HD 149404 


2012 Jul 04 


600 


C Oph 


2012 Jul 04 


176 


HD 152236 


2012 Jul 04 


240 


Elias 29 


2012 Jul 04 


3600 


AFGL 2136 IRS 1 


2012 Jul 06 


3960 


HD 179406 


2012 Jul 12 


1200 



Notes: The telluric standard star HR 6508 was 
observed on 2012 Aug 02 with an integration time 



of 336 s. 



TABLE 2 Absorption Line Parameters 









VLSR 


FWHM 


w x 


a(W x ) 


N(J) 


a(N) 


Target 


Molecule 


Transition 


(km s" 1 ) 


(km s" 1 ) 


(10~ 7 /xm) 


(10~ 7 urn) 


(10 13 cm" 2 ) 


(10 13 cm" 2 ) 


X Oph 


HF 


mo) 




3 




0.34 




0.31 


HD 149404 


HF 


R(0) 




3 




0.24 




0.22 


HD 152236 


HF 


R(0) 




3 




0.37 




0.34 


HD 154368 


HF 


R(0) 


4.9 


3.7 


1.80 


0.41 


1.66 


0.37 


HD 154368 


HF 


R(l) 








0.17 




0.23 


Elias 29 


HF 


fltoj 


3.0 


4.0 


19.5 


0.70 


18.0 


0.65 


Elias 29 


HF 


R(l) 








0.19 




0.27 


AFGL 2136 IRS 1 


HF 


B0) 


20.3 


11.6 


26.2 


1.64 


24.2 


1.52 


AFGL 2136 IRS 1 


HF 


R(l) 


24.5 


14.8 


17.8 


3.10 


24.9 


4.34 








fLSR 


FWHM 


Wx 


a(Wx) 




a(N) 


Target 


Molecule 


Transition 


(km s" 1 ) 


(km s" 1 ) 


(1Q- 7 nm) 


(10~ 7 (im) 


(10 16 cm" 2 ) 


(10 16 cm" 2 ) 


AFGL 2136 IRS 1 


H 2 


v \ 96,3^85,4 


23.7 


13.4 


30.0 


2.04 


3.38 


0.23 


AFGL 2136 IRS 1 


H 2 Q 


V\ 107,4-96,3 


24.8 


17.6 


21.1 


3.36 


2.43 


0.39 



Notes: Columns 4 and 5 are the line-center velocity and velocity full-width at half-maximum (including instrumental 



broadening effects) found by a Gaussian fit to the absorption feature. For non-detections of the -R(O) line a FWHM of 
3 km s _1 is adopted in determining the uncertainty in the equivalent width. For non-detections of the R(l) line, the 
observed FWHM of the i?(0) line is adopted in determining the uncertainty in the equivalent width. Columns 6 and 7 
are the equivalent width, W\ , and its la uncertainty. Columns 8 and 9 are the column density in the lower state, N, 
and its la uncertainty. HF column densities were de rived in the optical ly thin limit using transition dipole moments 
o f |M«(0)| 2 = 0.0104 D 2 and |/i fl( i)| 2 = 0.00693 D 2 HZemke et aLlUMlTl . and transition wavelengths of 2.499385 /an 
and 2.475876 fim (Webb & Rao 1968) for the R(0) and R(l) transitions, respectively. Note the difference in units for 
the HF and H 2 column densities. 



TABLE 3 HF and H 2 Abundances 

A r (HF) 7V(H 2 ) H 2 AT(CH) CH JV(HF)/JV(H 2 ) JV(HF)/JV(CH) 
Target (10 13 cm" 2 ) (10 21 cm' 2 ) Reference (10 13 cm' 2 ) Reference (10~ 8 ) 

X Oph < 0.93 0.43 ±0.10 1 3.02 ± 0.21 6 < 2.18 < 0.31 

HD 149404 < 1.99 0.62 ± 0.05 2 2.88 ±0.20 6 < 3.23 < 0.69 

HD 152236 < 2.02 0.53 ±0.13 2 2.40 ± 0.20 6 < 3.79 < 0.84 

HD 154368 1.66 ± 0.37 1.44 ± 0.40 3 6.31 ± 0.44 6 1.15 ±0.41 0.26 ± 0.06 

Elias 29 18.0 ± 0.65 26 ± 5 4 ... ... 0.69 ±0.14 

AFGL 2136 (LoS) 49.1 ±4.60 56 ± 9 4 ... ... 0.88 ±0.16 

AFGL 2136 (warm) 162 ±28.2 55-98 5 ... ... 1.65-2.94 

AFGL 2136 (cold) 15.0 ± 2.21 26-46 5 ... 0.33-0.58 

Notes: Upper limits on A (HI') in a single cloud component arc taken to be 3 times the la uncertainties reported 
in Table [2] Upper limits on the total line-of-sight column densities of HF (column 2) acc ount for mult iple clouds 
(velocity components) observed in CH (3 components for HD 149404, and 2 for HD 152236; Grcdcl ct al. 1993), and 
are used in calculating upper limits to the HF/H 2 and HF/CH abundance ratios. Upper limits on the J = 1 column 
densities toward HD 154368 and Elias 29 are small compared to observed J = column densities so we set the total 
HF column densities equal to the latter. For AFGL 2136 IRS 1 we present results for multiple cases described in this 
paper. The observed line-of-sight column densities (row labeled LoS) are presented, where the HF column density is 
taken to be t he sum of the columns in the J = and J = 1 levels and the H 2 column densit y is that reported by 
IKulesal 12002T ). Rows labeled " warm" and "cold" present the results from our analysis in Section l4. 21 where we assume 
both warm and cold components are present. This decomposition into warm and cold gas components is likely the 

more realistic pic ture. 

References: (l) ISavage eTaLI JT977T) ; (21 IRachford et al.l (I2009T) ; (3) IRachford et al.l ((2003 1 ; (4) IKulesal p002T ) (Note 
that only the para-H 2 column density was measured — via the v = 1—0, S(0) transition — and that the total H 2 column 
density was estimated assuming thermal equilibrium bet ween ortho and para H 2 , and a kinetic temperature derived 
from GO observations also presented therein); (5) Section 14.21 herein; (6) D. E. Welty (2002, private communication) 



